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In recent years localized waves have aroused great interest in the scientific community pj]. In relation to this topic, 
many efforts were devoted to the analysis of Bessel beams because of their unusual features: they are non-diffracting 
0, H 0, S @] and show superluminal behavior both in phase and group velocities 0, [1, 0, [13, [Hi • 

The aim of this paper is the study of the propagation of a Bessel beam through two absorbing layers, limited by 
two different half-spaces. Our approach will be based on the scalar analysis, since this analysis was proved to be an 
excellent approximation of the vectorial field which describes a Bessel beam [l2j . 

Let us consider a system formed by two layers, 1 and 2, limited by two half-spaces denoted by and 3, as sketched 
in Fig. [11 

All planes limiting the four different media are parallel and of infinite length. A localized wave, at frequency 
impinges on medium 1 coming from medium 0, which is non absorbing. On the contrary, media 1, 2, and 3 are 
absorbing. Be di and ^2 the thickness of the layer 1 (medium 1) and 2 (medium 2) respectively. Let us denote with 
the subscripts i, r and t the incident, reflected and transmitted fields, respectively. The subscripts 0, 1, 2, and 3 
characterized the media in which the field propagates. Be ng, ni, ?i2, the refractive indices of media 0, 1, 2, 3, 
respectively. For the time being, let us consider an ideal perfect homogeneous system. 

Let us consider a localized wave of Bessel type which impinges in medium 1 at normal incidence. We wish to 
remind that a Bessel wave - or Bessel beam - originates by the interference of an infinite number of plane waves, each 
of which tilted of the same angle 6q with respect to a given direction, say z, which is the direction of propagation 
of the beam. We refer to a specific Bessel beam, namely the one described in Ref. ^12]. For a Bessel beam like the 
one described in this particular reference, the field is linearly polarized and the scalar approximation is justified since 
two components of the field are negligible as compared to the third component, which is the dominant one. With 
reference to a cylindrical coordinate system (p, ipj z), a Bessel beam E is given by 
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where Jq is the Bessel function of first type, 6*0 is the axicone angle which characterizes the Bessel beam, n is the 
refractive index of the medium, and ka — uj/c is the wave number in vacuum (the field E is rotationally symmetric 
and thus independent on the angular coordinate Thus, the electric fields in media 0, 1, 2, 3 can be written as (see 
Fig. [U 
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FIG. 1: Schematic representation of propagation. A Bessel beam impinges on medium 1 at normal incidence, then in medium 
2 and, lastly, propagates through medium 3. Media 1,2 and 3 are absorbing media. 
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Eor = ro Joe-*'=°"°^'=°"^», (z < 0) (3) 

= tiJoie*'=°"i^™''^\ (0<z<di) (4) 

= n Joie-*'"i"™'^^\ (0<z<di) (5) 

£^2t = t2 Jo2 e*'=«"^(--'^l)™^»^ (di < z < (di + da) (6) 

i^2r = r2 Jo2 e-^'=«"^(^-'^l)=°''^^ (di < z < (di + da) (7) 

£;3 = t3Jo3e*'=«"^(^-(''i+''^»™"^% (z>(di+d2)) (8) 



where 

Jo = Jo(fco?^opsin6'o), Joi = Jo(fconipsin6'i), J02 = Jo(fcon2psin6'2) J03 = Jo(fcon3psin6'3), 

and 01, 62, and ^3 are the complex angles which characterize the Bessel function inside the media 1, 2 and 3, 
respectively. The temporal factor e~*"', which is present in all fields, is omitted for the sake of simplicity. 

All the refractive indices are complex, with the exception of no. The quantities r, t, which are labeled with reference 
to the medium in which the field propagates, denote the complex reflection and transmission coefficients, respectively. 

It is crucial to note that, even if the refractive index ni and the angle 9i are both complex, the product ni sin^i is 
real p^. In fact, in the propagation of a plane wave in a conducting medium, the surfaces of constant amplitude are 
parallel to the interface [14]. Thus, in order to meet this condition the product ni sin^i has to be real, otherwise the 
constant-amplitude surfaces would be plane, but with an inclination with respect to the interface. By considering the 
continuity of the phase of the incident and reflected waves, it is also possible to demonstrate that ni sinOi = no sin 

0: . . 

Since all the products n„ sin0„ are real and are equal to no sin^Oi we come to the remarkable conclusion that the 
shape of the Bessel beam is not modified when propagating into an absorbing medium. As a consequence, Eqs. (H])-® 
can be written as 



Eoi = j(,e''=""o^coseo^ (9) 

Eor = ro Joe-*'^'"""^™"^", (10) 

Eit = ii Joe*'^'""!^™"''!, (11) 

Eir = n Joe-*'=''"i^™^''\ (12) 

E2t - ^2 jQe*fco"2(^-<ii)cose2^ (^3^ 

E2r = r2 Joe-*'=''"^(^-'^l)™"^^ (14) 
£^3 = Joe*'=°"-^(^-(''i+''^»™''^^ (15) 

In order to evaluate the reflection and transmission coefhcients, we shall apply the boundary conditions to the field 
(that is the condition of continuity of the field on the border surface between two media), and to its first derivative 
with respect to the direction of propagation. The latter condition works only in the scalar approximation, and is none 
other than the condition of continuity of the tangential component of the magnetic field. In formula we have: 
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FIG. 2: Intensity of the reflected field in medium (Eq. [Sj, as a function of the radial coordinate p. Near p ~ the intensity 
of the Bessel beam is greater than that of the plane wave. On going away from the origin, the intensity of the beam decreases 
and the intensity of plane wave surmounts. Parameter values are: 6^0 = 30°, nir = 1-1, n2r = 2.2, d\ = 25 cm, di = \ cm, v = 2 
GHz (a), and v = ?> GHz (b). The values of the refractive indices, nir and nir, refer to dry sandy ground (medium 1) and glass 
(medium 2), respectively. The n2r = 2.2 value may also refer to other kinds of materials, such as stone or ceramics materials 



ai = fcorT-idi cos ^1, q;2 = fcp 77,2^2 cos 6*2, 

00 = fco'^o cos ^0? 01 = fepni cos 01, 02 = ^0*^2 cos 02j 03 = fco^-s cos ^3, (18) 
Eqs. (fT6|) and ((T7)) can be written as: 



1 + To = + ri 

iie'"i +rie-*"i =t2+r2 (19) 
^3, 



126 + r^e. — 13 



0o(l - ro) = 0i(ii - ri) 

01 (tie^i - rie-^"i) = 02(^2 - ^2) (20) 

02 (i2e^"^ -r2e-^"=) = i303. 

By solving equations (fTO|) and (|^D|) . we obtain the reflection coeflicient tq related to the propagation in the incoming 
medium 0, and the reflection and transmission coefficients related to the propagation in the layers 1,2 and in the 
half-space 3 (the calculations are rather laborious but present no difficulty) 

a0o-ai0i 
acfiQ + ai (pi 



^ ^ 2016-'°^ (l + rp) 
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and 

^+ = 1 + 6-2-2$, A- = 1-6-2-2$. 

By replacing Eqs. (|2i p - ((26|) in Eqs. ([9|)- (fT5|) . we are now able to obtain the electromagnetic fields that describe the 
propagation of a Bessel beam through the multilayer system of Fig. [TJ 



We are mainly interested in evaluating the reflected field in medium 0, in a physical situation in which medium 3 
is equal to the medium 1, that is for ni — n^. This situation is related to the possibility of using Bessel beams in 
searching for objects buried underground as well as any kind of macro impurity or defect inside a given material. 

Because of the interference of multiple reflections on the layers, the amplitude of the reflected field is an oscillating 
function with respect to the frequency. For this reason, the intensity of the beam greatly depends on the frequency. 
Thus, for given parameter values, the choice of the frequency is extremely important in order to have the Bessel beam 
amplitude surmounting the amplitude of the plane wave. 

In Fig. [2] we show the intensity of the refiected field (Eq. [3]) from a Bessel beam and a plane wave as a function of 
the radial coordinate p, for two different frequencies. The figure is limited to the central portion of the Bessel beam, 
since it is the only one of physical interest. 

In order to check the validity of the scalar approximation, for the value of the parameter 9o used here {9q = 30°), 
and for p in the range of physical interest, we have verified that one component of the electric field is much larger 
than the other. The components of the field became comparable only close to the first zero of the Bessel function, 
and, in this case, the scalar approximation loses its validity. 

The fields which describe the propagation of the plane wave through the layers are given by Eqs. ((21)-® by putting 
9„ = 0. 

Our parameter values refer to the case of a sheet of glass buried underground |15j]. Near p — the gain is evident 
for both frequencies. Even if the intensity of the reflected field is only a few percent with respect to the incoming field 
[l^ . the intensity of the Bessel beam in both cases is about three times greater than that of the plane wave. 

For complex refractive indices, we chose a mean standard value in the microwaves range. Moreover, we set the 
value of the imaginary part of the dielectric constant at approximately one order of magnitude less than the real one. 
We are aware that this position may appears to be rather rough. However, it is a plausible approximation for our 
purposes, since our aim is to analyze the difference in the propagation between a Bessel beam and a plane wave, 
rather than to analyze a specific material (under specific conditions of humidity, temperature etc.). We should recall 
that the the complex refractive index Uc = nr + iui is related to the complex dielectric constant Cc = e + ie by means 
of the relation — e^- Thus, for e = 0.1 e we simply have rii ~ 0.05 n^. 

Once the use of Bessel beams has been shown to be advantageous with respect to that of plane waves, we can 
investigate the difference of intensity in the refiected field in the presence or absence of a given buried material. 

In Fig. [3] we report the three-dimensional intensity of the reflected field as a function of the Cartesian coordinates x 
and y {p = y^x^ + y^). The higher signal refers to the passage in the presence of glass, while the lower one refers to the 
passage only through the ground. The difference in the intensities is evident so that, in an experimental investigation, 
the presence of a buried object (glass, stone or ceramics materials) should be clearly detected. 



We have shown that the use of Bessel beams can be advantageous as compared to plane waves, in the detection of 
buried objects. Our analysis refers to a particularly favorable case, that is, objects buried inside very dry material 
(e.g. sandy or clayey ground) |,17i |. In this case, in fact, microwaves undergo small absorption. 

We showed that, near p = 0, Bessel beams always present a gain with respect to plane waves, provided that 
a suitable frequency is chosen: the gain becomes more and more negligible as the frequency increases. For high 
frequency values, the advantage of using a Bessel beam is lost, since the amplitude of the plane wave always exceeds 
that of the beam. 
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FIG. 3: Three-dimensional intensity of the reflected field (Eq. |3]) as a function of the Cartesian coordinates x and y. Higher 
signal refers to the reflected field in the presence of buried object; smaller signal (dark gray) refers to the propagation through 
the only ground. Parameter values are: do = 30°, di — 25 cm, d2 = 1 cm, and (a): nir = 1.2, n2r = 2.25, v = 2.2 GHz; (b): 
nir = 1.4, n2r = 2.3, u = 2.32. We now need to make a comment or two on the size of the Bessel beam generator: for 6o = 30°, 
in order to have a field depth of 50 cm, the converging system must have have a radius of about 30 cm This dimension is 
not unusual for an electromagnetic mirror, while a lens of this size would have to be custom-made. Naturally, the size of the 
microwave generator-launcher is different, depending on the frequency. 



A second, and perhaps more important, aspect related to Bessel beams is that the difference between the propagation 
through the ground alone and the propagation through the ground in the presence of buried objects is quite evident 
(Fig.© [p. 

In the light of these conclusions, we can think that a GPR (ground penetrating radar) system which utilizes a 
Bessel beam as incoming pulse could present some advantages as compared to a traditional GPR apparatus. 

Our theoretical model works under two main approximations: namely, the smooth surfaces of the layers and 
homogeneous media. From an experimental point of view, the first approximation works well as long as the dimensions 
of roughness are sufficiently small as compared with the wavelength; the second one requires a specific software, in 
order to select between signal and noise. 

Because of their feature of localized waves (i.e. localized energy), Bessel beams could also provide further information 
about the dimension and shape of buried objects. However, an analysis devoted to this aspect goes beyond the purpose 
of the present work and will be reported elsewhere. 
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